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Different methods have been recently put forward and implemented experimentally to inverse en- 
gineer the time dependent Hamiltonian of a quantum system and accelerate slow adiabatic processes 
via non-adiabatic shortcuts. In the "transitionless tracking algorithm" proposed by Berry, short- 
cut Hamiltonians are designed so that the system follows exactly, in an arbitrarily short time, the 
approximate adiabatic path defined by a reference Ifamiltonian. A different approach is based on 
designing first a Lewis-Riesenfeld invariant to carry the eigenstates of a Hamiltonian from specified 
initial to final configurations, again in an arbitrary time, and then constructing from the invari- 
ant the transient Hamiltonian connecting these boundary configurations. We show that the two 
approaches, apparently quite different in form and so far in results, are in fact strongly related 
and potentially equivalent, so that the inverse-engineering operations in one of them can be rein- 
terpreted and understood in terms of the concepts and operations of the other one. We study as 
explicit examples the expansions of time-dependent harmonic traps and state preparation of two 
level systems. 

PACS numbers: 37.10.De, 32.80.Qk, 42.50.-p, 03.65.Ca 



I. INTRODUCTION 

Externally imposed time-dependent interactions are 
frequently varied slowly to keep adiabaticity and control 
the final state of a quantum system robustly versus pa- 
rameter fluctuations. There are however many instances 
where we would like or need to quicken these operations. 
If they take too long, they are impractical for applica- 
tions in which they are repeated many times, e.g. to 
manipulate quantum information and transport ions or 
atoms or they may suffer from decoherence, noise 
or losses, so speeding them up may be the only way to 
actually implement the ideal final outcome. Moreover in 
many experiments, as in atomic fountain clocks, high rep- 
etition rates contribute to achieve better signal-to-noise 
ratios and better accuracy 0. Adiabatic steps are also 
the bottleneck in some cyclic processes, and determine, 
for example, the cooling rates in quantum refrigerators 
quantifying the unattainability of absolute zero [j-Q- 

Recently, several works have been devoted to theoret- 
ical proposals 0, or experimental realizations 
[22j of fast non-adiabatic shortcuts to the states reached 
by a slow adiabatic processes in matter wave expan- 
sions or compressions, splitting, and transport. Berry, 
in particular, has proposed a "transitionless tracking al- 
gorithm" to design time dependent interactions so that 
the system follows exactly, in an arbitrarily short time, 
the approximate adiabatic path defined by a reference, 
zeroth order Hamiltonian [l^]- This method has been 
applied to speed up adiabatic passage techniques and 
achieve fast and robust population control in two and 
three level atomic systems jl^]. A different approach is 
based on designing first a Lewis-Riesenfeld invariant jlj 
to carry the eigenstates of a Hamiltonian from initial 
to final configurations, again in an arbitrary time, and 
then constructing from the invariant a transient, driving 
Hamiltonian 0, [I]- These methods were compared for 



harmonic oscillator expansions for which they provided 
rather different shortcut paths 

Berry mentioned the existence of connections between 
the invariants and the transitionless algorithm for a two- 
state system without pursuing them further [loj . Fol- 
lowing that hint we show in this work that these two 
approaches are in fact closely related, and can be stated 
in common terms, so that the inverse-engineering opera- 
tions in one of them can be reinterpreted using concepts 
and operations of the other one. This sets their poten- 
tial equivalence. The different results that have been 
found so far, as explained in detail below, are due to the 
ample freedom offered by both approaches to construct 
the driving shortcut Hamiltonian. We study the general 
setting as well as two explicit examples: expansions of 
time-dependent harmonic traps and state preparation in 
two level systems. 

A word of caution on notation: as we shall deal with 
different methods and examples, multiple usage of some 
symbols such as l^n), I"-}, H, I, A„, U , or En is unavoid- 
able unless we load them with sub and superscripts, so 
consistency is strictly guaranteed only within each sec- 
tion. In most cases this repeated usage will suggest a 
possible relation. The context and explanations will clar- 
ify how this comes about. 



II. GENERAL FRAMEWORK 
A. Lewis-Riesenfeld invariants 

We shall describe first Lewis-Riesenfeld theory in a 
nutshell [24| . Let us consider a quantum system evolving 
with a time-dependent Hamiltonian H{t). A dynamical 
invariant I{t) satisfies 



(1) 
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so that its expectation values remain constant in time. 
I{t) can be used to express an arbitrary solution of the 
time-dependent Schrodinger equation 



(2) 



as a superposition of "dynamical modes" |'0n(^))i 

= ^c„|^„(t)), (3) 

n 

\Mt)) = e^""(*)|0„W), (4) 

where c„ are time- independent amplitudes, \(f>n{t)) is 
the n-th eigenvector of the invariant I{t), = 
An|0n(i)), with A„ constant, and the Lewis- Riesenfeld 
phase is defined as [l^] 



ih 



d_ 



H{t') Mt'))dt'. (5) 



We use for simplicity a notation for a discrete spectrum 
of I(t) but the generalization to a continuum or mixed 
spectrum is straightforward. We also assume a non- 
degenerate spectrum. 



B. Invariant based inverse engineering 

Suppose that we want to drive the system from an ini- 
tial Hamiltonian H{0), to a final one H{tf), such that the 
populations in the initial and final instantaneous bases 
are the same, but admitting transitions at intermediate 
times. To inverse engineer a time-dependent Hamilto- 
nian H(t) and achieve this goal, we may define first the 
invariant through its eigenvalues and eigenvectors as 



iit) = J2\Mt))K{Mt)\ 



(6) 



The Lewis-Riesenfeld phases may also be chosen as arbi- 
trary functions to write down the time-dependent unitary 
evolution operator [/, see Eqs. pi4p . as 



It must obey 



ih^U = H{t)U, 



(7) 



(8) 



which we solve formally for H(t), 

H{t)^ih{dtU)U\ (9) 

that is. 

Hit) = F{t) + ihY, \dt(jyn{t)){Mt)\, (10) 



where Fit) is diagonal in the basis of the invariant, 

F{t) - -ri^|0„(t))d„(0„(t)|. 



(11) 



and the dot denotes derivative with respect to time. Note 
that for a given invariant there are many possible Hamil- 
tonians corresponding to different choices of phase func- 
tions an(i).l23| In general /(O) does not commute with 
i/(0), which means that the eigenstates of /(O), |0„(O)), 
do not coincide with the eigenestates of (0). Hitf) does 
not necessarily commute with I{tf) either. If we impose 
[/(0),7?(0)] = and [I{t}),H{tf)] = 0, the eigenstates 
coincide and then a state transfer without final excita- 
tions is guaranteed. In a typical application the Hamil- 
tonians H(0) and H{tf) are given, they set the initial 
and final boundaries for the process, and we use these 
boundary conditions to define lit) and its eigenvectors 
accordingly. A convenient, although by no means neces- 
sary, relation is to set /(O) = H{0). We shall see specific 
examples of how this works in Sees. Ill and IV. 



C. Transitionless tracking algorithm 

In Berry's method [lol |. the starting point is a time- 
dependent reference Hamiltonian 

i^oW = El"oW)^"°'(t)(no(t)| (12) 

n 

for which the approximate time-dependent adiabatic so- 
lutions are 



|V'„(<))=e'«"(*)|no(t)) 



(13) 



where the adiabatic phases, with dynamical and geomet- 
ric parts, are 

Ut) = -l f dt'El^\t') + i f dt'{noit')\dt'no{t')). 

" Jo JQ 

(14) 

The approximate adiabatic vectors in Eq. (|13p are de- 
fined differently from the dynamical modes of the pre- 
vious subsection, but they may potentially coincide, as 
discussed below, so, with some caution, we use the same 
notation. Defining now the unitary operator 



U = Ee^«"(*Vo(t))(no(0)|, 



(15) 



a Hamiltonian Hit) can be constructed, using again the 
general relation to drive the system exactly along the 
adiabatic paths of Ho{t) as 

H(t) = Hoit)+Hiit), 
Hi{t) = ihY, (\dtno{t)){no{t)\ 

n ^ 

- {no{t)\dtno{t))\no{t)){no{t)\) , (16) 
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where Hi{t) is purely non-diagonal in the {|fio(i))} basis. 

We may change the functions En\t), responsible for 
the dynamical part of the phase, and therefore Ho{t) it- 
self, keeping the same \no{t)) eigenvectors. We could for 

example make all the En\t) zero to suppress the dy- 
namical phases, or compensate the geometric phase to 
have = Therefore, the Hamiltonian can be 

generally written in terms of the phases as 

H{t) = G{t) +inY^ \dtno{t)){no{t)\, (17) 



where 



G{t) = -hJ2\noim„(noit)\ 



(18) 



is diagonal in the instantaneous basis of Hq. Subtracting 
Hi, Hq may also be written as 



In general it is required that Hi{t) vanish for t < and 
i/ > 0, either suddenly or continuously at the extreme 
times. In that case the \no{t)) become also at the extreme 
times (at least at t = 0~ and t = f^) eigenstates of the 
full Hamiltonian. 

Using Eq. ([ij and the orthonormality of the {|no(0))} 
it is easy to check that we may write invariants of II[t) 
with the form 



/(i) = ^|rio(i))A„(no(i)l, 



(20) 



where the A„ are constant eigenvalues. For the simple 

choice A„ = E^n\o), then /(O) = i?o(0). 

Up to now we have presented the invariant-based 
and tracking algorithm approaches in a common man- 
ner to make their relations obvious. By reinterpreting 
the phases of Berry's method as ^„(t) = a„(i), and the 
states as |no(i)) = \(j)n{t)), we may immediately equate 
G{t) = F{t) and the Hamiltonians H{t) in Eqs. ^ 
and ([T7|) . We may also find the Ho{t) implicit in the 
invariant's method using Eq. p9p . In other words, the 
dynamical modes in the invariant-based method can be 
also understood as approximate adiabatic modes of a cer- 
tain Hamiltonian Ho{t). 

An important caveat is that, although the two meth- 
ods could coincide, they do not have to. Given H{0) 
and H(tf) there is much freedom to connect them us- 
ing different invariants, phase functions, and reference 
Hamiltonians Hq^I). It should be clear by now that each 
of these methods does not provide a unique shortcut but 
entire families of them, a welcome flexibility that allows 
to optimize the path according to physical criteria and/or 
or operational constraints. 

In the following sections, we shall work out two specific 
examples where the connections, differences and similar- 
ities of the two approaches are illustrated and examined 
further. 



III. TIME-DEPENDENT HARMONIC 
OSCILLATOR 

A. Lewis-Riesenfeld invariants 

The Hamiltonian of a particle in a time-dependent har- 
monic oscillator with varying angular frequency (as all 
"frequencies" hereafter) w(i), is 



H{t) 



2m 



P" 



muj-^ (t) 



(21) 



The instantaneous eigenstates and energies are, respec- 
tively. 



{x\n{t)) 



f^W 1/4 

mu!(t) \ 



1 



(2»n!)i/2 



X exp 



muj{t) 2 



muj{t) 



,(22) 



where if„ is a Hermite polynomial, and En{t) = (n + 
l/2)hu{t), so that H{t)\n{t)) = EJt)W{t)). Trying a 
quadratic ansatz for the invariant |24l427l |. 



m 



TO 



(23) 



where tt = b{t)'p — mbq plays the role of a momentum 
conjugate to q/b, and inserting it into Eq. ([1]), the scaling 
factor b = h{t) is found to satisfy the Ermakov equation 



6 + a;2(t)&: 



Wo 
63 



(24) 



is in principle an arbitrary constant, which we fix as 
the initial frequency. The eigenfunctions of I{t) are 



{x\4>n{t)) 



1 



(2"n!6)i/2 
xH, 



exp 



2h \ b 



\~ir) b 



(25) 



and, since I{t) has the structure of a generalized har- 
monic oscillator, i.e., quadratic but with "crossed" 
momentum-position terms, the eigenvalues are A„ = 
{n + l/2)huJo. 

Substituting the wave function {x\(l)n{t)) into Eq. (jlOp . 
the Hamiltonian may be written as 



b , , TO /-o UJ, 

H{t) = m + ^{pq + 9P) - y ( & + ) 9 



2 ^ 

\ ~2 



(26) 



For consistency with Eq. (j2ip the crossed terms must 
cancel. This is indeed the case. F{t) can be determined, 
using Eq. (IllL from the phase of the dynamical modes, 
Eq. ([5]), which is now 



(27) 



4 



This gives d„ = ~{n + \)uJo/b'^ and F{t) = J(t)/6^ 
Hamiltonian H{t) can be finally written as 



H{t) 




The 



(28) 



which is nothing but the Hamiltonian (I^Tl) after substi- 
tuting using the Ermakov equation (p4|) . 



B. Invariant- based engineering approach 

In the inverse engineering approach based on Lewis- 
Riesefeld invariant theory as presented in the main 
goal is to find a "trajectory" for the external parameter 
Lij{t) SO that the populations of the final oscillator levels 
are the same as the populations of the initial one. De- 
signing I{t) first here means to design b{t) . Let us assume 
an expansion with initial and final frequencies w(0) = ojq 
and i^{tf) — LOf. To make I{t) and H{t) commute at 
t = and t f and have common eigenfunctions we impose 
the boundary conditions 

6(0) = 1, 6(0) = 0, 6(0) = 0, 



ojf/ujo, bitf) = 0, b{tf) - 0, (29) 



by comparing Eqs. (p2|) and (l25|) and using Eq. (1241) . 
A simple polynomial ansatz can be used to interpolate 
h{t) at intermediate times [1]. Once b{t), and therefore 
the invariant, are set, the time-dependent frequency uj[t) 
and the shortcut Hamiltonian follow from the Ermakov 
equation (1^ . This method, including the effect of grav- 
ity, has been realized experimentally 2l|, and extended 



to Bose-Einstein condensates 0, H^l , and to design trans- 
port protocols jl7l] . 



Because of the boundary conditions (|29l) , Hi it) vanishes 
at i = and t = tf. The crossed terms are canceled out 
in the full Hamiltonian H(t) = Ho{t) + Hi{t) given by 
Eq. 



D. Transitionless tracking algorithm (standard 
application) 

Unlike the previous subsection, in a more standard ap- 
plication of the transitionless tracking algorithm to the 
harmonic oscillator [ll|, Ho{t) is set first as an ordi- 
nary harmonic oscillator with given frequency uj(t), i.e., 
it takes the form so that \n{t)) in Eq. should be 
now reinterpreted as \nQ{t)). Hi(t) is calculated from Eq. 
(I16L and the resulting shortcut Hamiltonian becomes 



H{i) = ^ + \mLo\t)e --^^ 
^ ' 2m 2 ^ ^ 4w(t) 



(pq + qp), (32) 



Ho{t) 



Hi{t) 



a generalized harmonic oscillator with crossed terms that 
imply a non-local interaction [26l |. 



E. Relation to invariants 

To reinterpret the previous subsection in terms of 
invariants, we construct I{t) = X^n l"o(i))'^n(?^o(^)|, 
with A„ = eI?{Q) = (n + l/2)fta;o- Since Ho{t) = 
J2n\''^oit))En \t){nQ{t)\ with instantaneous eigenvalues 

E^i\t) = [n + l/2)hLu{t), the invariant is now propor- 
tional to an ordinary harmonic oscillator. 



Wo 



lit) = -^Mt) 



^0 



2^ + 2"^" 



(33) 



C. Invariant-based method in 
transitionless-tracking-algorithm language 



As pointed out in Sec. IHCi the invariant based method 
can be restated in the language of the transitionless track- 
ing algorithm. To find the reference Hamiltonian 7Jo(i) 
implicit in Sec. IHI B[ we interpret the |'/)n(i)) as the 
eigenstates |f^o(^)) of Hq, take ^„(i) = ck„(t), and set the 
eigenvalues according to Eq. p^ . This gives 




bb 



2^0 



(30) 



which is a generalized harmonic oscillator. The corre- 
sponding Hamiltonian Hi(t) takes the form 



rr . ^ b m 9 6M . 



b^-bb 



(31) 



and /(O) = i/o(0). Using = -Ei°\t)/h + 
i{nQ(t)\dtnQ{t)), and {no{t)\dtnQ{t)) = 0, then letting 
(^nit) = init) and \4>nit)) — \no(t)), we may write down 
H{t) from Eq. (HU]), 



Hit) = Fit) - 



Aioit) 



ipq + qp), 



(34) 



but here Fit) = Hoit) — [a;(t)/ajo]-f (i); so we recover the 
Hamiltonian /(i) does not commute with Hit) in 

general. To guarantee that /(t) and Hit) have common 
eigenstates at i = and t = tf, the boundary conditions 
w(0) — and w(t/) = should be satisfied, so that i?i(t) 
vanishes at the initial and final times, as discussed in Sec. 

El 



IV. TWO-LEVEL ATOM 

The two-level atom is another fundamental model. 
Speeded-up adiabatic state preparation methods such as 



5 



Rapid Adiabatic Passage (RAP) in two-level atomic sys- 
tem may be useful in chemical reaction dynamics, laser 
cooling, or quantum information processing. 

For the two-level atom, using |1) = (°), |2) = (g), 
the time-dependent Hamiltonian which we consider, in a 
laser adapted interaction picture and applying the rotat- 
ing wave approximation, is 



(35) 



where A = A{t) and flu = il.R{t) are the time-dependent 
detuning and Rabi frequency, and cp = ip{t) a time- 
dependent phase. The instantaneous eigenvectors are 



\n+{t)) = cos( -je'n2)+sin( -)|1), (36) 



\n.{t)) = sin(^-j|2)-cos(^-je-^'^|l), (37) 

with the mixing angle 9 = 9{t) = arccos(A/ri) and eigen- 
values E±(t) ^ ±hn{t)/2, where n{t) = ^A^ + n%. If 
if = 0, and the adiabaticity condition 



QrA - fiffA 



173 



«1, 



(38) 



is satisfied, the state evolving from |V'±(0)) — \n±{0)) 
follows the adiabatic approximation 

l^±(i)) =exp|-^ f dt'E±{t')\\n±{t)), (39) 



whereas transitions will occur otherwise. If the state 
starts from |n+(0)), the adiabatic evolution of the popu- 
lation of levels 1 and 2 is 



Pf{t) = |(lK(t))P = sin2 - , 



= |(2K(t))|^=cos^ - 



(40) 



In what follows we shall speed up the adiabatic passage 
by the two methods presented here. 



A. Invariants method 

To use the invariant-based inverse engineering method, 
we first parameterize the eigenvalues and eigenstates of 
the invariant /(i), satisfying /(i)|(/)„(i)) = A„|0„(i)) con- 
sistently with orthogonality and normalization, in paral- 
lel to Eqs. ((361) and (jST)) . 



\Mt)) 
\^-{t)) 



;(|)e^^|2)+sin(|)|l), (41) 
|2)-cos(^)e-'^|l), (42) 



and set \± ~ ±.MIq/2. Thus, I{t) can be expressed as 



ut) ^ -no 

^ ^ 2 V sm7e 



cos 7 sm7e 

1/3 



1/3 

cos 7 



(43) 



where /3 — /3{t) and 7 = "f{t) are auxiliary time- 
dependent angles. Using Eqs. ([55]) and (HTI). the Lewis- 
Riesenfeld phase (O is now calculated as [2^ 



1 



±{t) = ±- / [A{t')-2n{t')]dt', 



which gives 



1 



a± 



A{t) - 2n{t) 



(44) 



(45) 



-I- sin 7 cos {/3 — (f). 



where 

n = {A + p)cos'{l) 
Substituting Eq. (gS]) into Eq. (HU]), we obtain 

'2n - a' 



Hit) = 



I{t)+ihJ2\dtMt)){<P±it)l (46) 



which can be finally expressed as 

M 



-M 



(47) 



where 



M — A cos^ 7 -t- sin 7 cos 7 cos {13 — (p) — $ sin^ 7, 
= [A cos 7 + sin 7 cos {(3 ~ ip) + $ cos 7] sin 7 — 17. 

The Hamiltonian (H71) must be equivalent to the Hamil- 
tonian dSSl), so from M = A and Ne'l^ = flRe"^ [H, we 
get the following auxiliary equations. 



7 = QRsm{l3 - ip), 



(48) 



(A -I- /3) sin 7 = n R cos J cos (13 - (p). (49) 

In general H{t) does not commute with I{t), 

[H{t),I{t)]/{n^nQ) = a+{Asm'ye'^ -nRCOSje'"^) /2 

- (7_ (Asin7e~''^ - 17ijcos7e~''^) /2 
-I- aoifl R sin J sin [13 — ip). (50) 

[iJ(0),/(0)] = is satisfied, if 

A(0) sin7(0)e*'3(°) - 17^(0) cos7(0)e''^(°) = 0, (51) 
A(0) sin7(0)e^*'3(°) - 17^,(0) cos7(0)e^"^(°) = 0, (52) 
nR{0) sin 7(0) sin [/?(0) - (^(0)] = 0, (53) 

and there are similar equations for tf. For population 
inversion processes we are interested in processes starting 



6 



and ending with zero flfi and some finite detuning, so we 
impose 



r!fl(0) = 0, 7(0) =n^, 



(54) 



so far with arbitrary /3(0) and (p{0). In this case, H{0) 
and 1(0) have common eigenvectors, which are exactly 
the pure ground state |1) and the excited state |2). 
Similarly, for [H {t f) , I {t j)] = 0, 



(55) 



with arbitrary f3{tf) and f{tf). Again H{tf) and I{tf) 
share common eigenstates, which are the ground and the 
excited state. 

Substituting the above boundary conditions into Eqs. 
and (P^ . we further obtain 



7(0) = 0, 7(t/) = 0, 



(56) 



whereas /3(0) and ${tf) will determine the value of the 
initial and final detunings. 

We are now ready to set some ansatz for 13 and 7 using 
appropriate boundary conditions. Once the functions (3 
and 7 are fixed, we can construct Qfj and A and thus 
the Hamiltonian H{t) with a given (p. In the following 
subsection, we will give some examples to show how the 
invariant-based engineering method works with different 
boundary conditions. 



B. Example 
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We shall apply the previous results to design fast popu- 
lation transfer protocols. For simplicity, we assume (p = 
and consider the, yet unknown, Hamiltonian H{t), 

where ^Iji and A are determined, from (|48| and (|49|) by 

Qr = j/sm/3, (58) 
A = cot7Cos/3 - /3. (59) 

We suppose that this Hamiltonian drives the state from 
|1) to |2), up to the phase factor, along the invariant 
eigenvector \(j)+{t)). To this end, we set the boundary 
conditions 17(0) = f2(f/) = and 

7(0) = TT, 7(0) = 0, (60) 
7(t/)=0, 7(t/) = 0. (61) 

As mentioned before, we have freedom to choose the val- 
ues of ^(0) and (3{tf). According to Eq. (1551) . it is useful 
to keep /3 close to (n-|- l/2)7r, so as to minimize ftR along 
the path, whereas the derivatives fix the initial and final 
detunings, see Eq. ([59|. which should have here opposite 
signs. Moreover they should not be too large to keep (3 



FIG. 1: (Color online), (a) Examples of polynomial ansatzs 
for -f{t) = Ej=oaji^ (solid red line) and /?(*) = E^^o^j^^ 
(dashed blue line), (b) The corresponding functions of i}u 
(solid red line) and A (dashed blue line) determined by Eqs. 
(|58|) and (|59|l . (c) Time evolution of the populations of lev- 
els 1 and 2: Pi (dashed blue line), P2 (solid red line), and 
adiabatic approximations P^'' (dotted purple line) and 
(dash-dotted orange line), see Eq. H40|l . hardly distinguish- 
able from the former, tf = Ifis. 



close to the chosen reference /3 value, and to avoid (3 — 
at some intermediate time and thus an infinite fiij. Con- 
sidering all these physical constraints, we impose 

/3(0) = -7r/2, /3(0) = 37r/(2i/), (62) 
/3{tf) = -n/2, /3(t/) = -37r/(2i;), (63) 

where the negative sign of (3, see Eq. ((55)1 and Fig. [I](a), 
keeps flu positive, as 7 becomes negative. 

To interpolate at intermediate times we assume a poly- 
nomial ansatz. Fig. [T] (a) shows j{t) — X]^=o ^i^'' 

f3{t) = J2^j=o ^3^'' ' where the coefficients are obtained by 
solving the equations set by the boundary conditions. 
The time-dependent ilji and A calculated from Eqs. ((58)) 
and ((55)1 are shown in Fig. [Tj (b) . Once we have speci- 
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FIG. 2: (Color online), (a) Examples of polynomial ansatzs 



for 7(f) = E, 



=0 '^j^'' 



(solid red line) and P{t) — E, 



(dashed blue line), (b) The corresponding functions of Qr 
(solid red line) and A (dashed blue line) determined by Eqs. 
(|58|l and (|59p . (c) Time evolution of the populations Pi 
(dashed blue line) and P2 (solid red line); adiabatic approx- 
imations Pi""* (dotted purple line) and P2'"' (dash-dotted or- 
ange line) for comparison, tf = Ips. 



f3{t) = ^^^QbjP. Fig. [5] shows the results, to be com- 
pared with those of Fig. [TJ Note that in Fig. [2] (b) Rabi 
frequency and detunings are smaller than in Fig. [T] (b), 
so smaller energies are involved. Now the dynamical evo- 
lution is not adiabatic, see Fig. [I](c). The method can be 
further complemented by optimizing the trajectory with 
resp ect to different physical cost functions or constraints 
|15| . for example, setting bounds for ilji, whose square is 
proportional to the laser intensity. This will be discussed 
elsewhere. 



C. Invariant-based method in transitionless 
tracking algorithm language 

In order to reexamine the invariant-based inverse engi- 
neering approach in the language of Berry's transtionless 
tracking algorithm, we take \(j)±{t)) as \no±{t)), and let 
a±{t) = C±(i), so that Eq. gives 

jj^^^^ ^ A cos 7 -t- flR sin 7 cos (/3 - ip) ^^^^ 



2n- A- 2/3cos2 (7/2) 



(66) 



and the Hamiltonian Hi (t) in Eq. (|16p is 



-/3sin^7 (-17 -I- I sin 27) 



{ij + 4 sin 27)6" 



13 sin 7 



(67) 

Using the boundary conditions ([5ll) - (|55)) . Hi{t) vanishes 
at i = and t = tf. 



D. Transitionless tracking algorithm 

Let us now apply Berry's transitionless tracking algo- 
rithm taking the Hamiltonian psp as reference Hamilto- 
nian Ho{t) [ill, 



fied H{t) in ((57)) . we solve the dynamics numerically by 
a Runge-Kutta method with adaptive step, see the pop- 
ulation inversion in Fig. [T](c) for levels 1 and 2. We have 
also compared the bare state populations Pi and P2 with 
the populations of the instantaneous eigenstates of H{t), 
Pf' and P^"*. Their agreement shows that the designed 
protocol is in fact an adiabatic passage for the specified 
final time tf. 

We may impose additional conditions at an intermedi- 
ate time, for example, to keep j3 closer to —tt/2, 

/3(0) = -7r/2, P{tf) = -7r/2, (3(tf/2)^-TT/2, (64) 
m = ^l{2tf), P{tf) = -iTl{2tf), (65) 

where we have also diminished the detuning. This new 
set of conditions requires a higher other polynomial. 



The driving Hamiltonian ([T5| becomes in this case 



Hi{t) 



-ifnin^e (-i0+ |sin26')e 

If SIT? 9 



2 \ {ie + sm2e)e-^^ 
For iy9 = 0, this reduces to 

^'^'^-2(-^na 



(69) 



(70) 



where = ^ = (fij^.A — 17;? A)/r2^. In [12|], this was used 
to speed up an Allen- Eberly scheme for Ho{t) and achieve 
fast population transfer. To see the physical meaning 
and realizability of the method we must go back to the 



8 



Schrodinger picture: For H{t) = Ho{t) + Hi{t) this im- 
plies using two lasers with the same frequency, orthogonal 
polarization, and time-dependent intensities but different 
intensity shapes [l2]. The alternative is to drive the sys- 
tem with Hi{t) only, without Ho{t). In the Schrodinger 
picture this amounts to act with one laser and to per- 
form level shift engineering to modulate the transition 
frequency so as to leave A = in the interaction picture 
[l^ . Note that these complications (an extra laser or the 
need for level-shift engineering) do not arise in the results 
obtained in the previous subsection. 

E. Relation to invariants 

Let us now reinterpret the previous (standard) Berry's 
transitionless tracking algorithm in terms of invari- 
ant theory. In the language of Lewis-Riesefeld invari- 
ant theory, we can construct an invariant as I{t) = 

E± \no±{t))\±{no±it)\, where A± = E^^\o) = ±hno/2, 
with matrix form 

^^^)-2^^\^i^ee-^^ -cose j' (^^^ 

and /(O) = HoiO). Since Ho(t) 

X]± l'^o±(i))^'±'(i)(^o±(OI with instantaneous eigenval- 
ues Ef\t) = ±hn{t)/2, we have i?o(0 = [fl{t)/no]I{t). 
Using e± = -E^l\t)/h + i{no±{t)\dtno±{t)) and 
(rio±(t)|9tno±(t)) = ±i(/5 cos^ (0/2), then letting 
l0±(i)) = Wo±{t)) sind a±{t) = S,±(t), we may write 
H{t) from Eq. ([TU|) . Canceling terms this gives exactly 
the Hamiltonian i?o(i) + Hiit) in the previous subsec- 
tion. I{t) does not commute with H{t) in general, but, 
when the boundary conditions 

0(0) = 7r, 0(0) =0, (72) 
e{tf)=Q, 9{tf)^0, (73) 

are satisfied, Hi(t) will vanish at initial and final times. 

V. CONCLUSION 

In previous publications we applied and compared two 
methods to speed-up adiabatic processes through non- 



adiabatic shortcuts: Berry's transitionless tracking al- 
gorithm and the invariant- based inverse engineering ap- 
proach. Their differences were emphasized, in particular 
in time-dependent harmonic oscillators or transport 
of particles by a moving trap [l3l- The message here 
is quite different, even opposite: we point out now that 
in fact both approaches share a common ground of con- 
cepts and structure. There is, however, no contradic- 
tion. It is indeed possible to interpret a particular in- 
verse engineering operation using either the language of 
transitionless-tracking or invariants approaches, and con- 
sider them to be potentially equivalent. The explanation 
of the differences found is the large freedom to design dif- 
ferent Hamiltonians for a given speed-up goal. In other 
words, the different results are not fundamental but due 
to the particular choices that have been made to resolve 
that freedom in specific implementations. The choice of 
method from this point becomes thus, in part, a mat- 
ter of taste, but there are also elements that make one 
or the other approach more natural or convenient. For 
example, systems with Hamiltonians that admit known 
structures for the invariants are easy to approach with 
the invariant-based method. This includes transport, ex- 
pansions, rotations psj . or, as shown here, discrete level 
systems. The tracking algorithm can be applied in many 
systems where the invariants are unknown. In summary 
this work provides a significant step towards a deeper 
understanding of shortcut-to-adiabaticity methods that 
will help to choose the most adequate approach in atomic 
transport, quantum gates, and generally atomic manip- 
ulation and control applications. 
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